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Abstract 

Analytical  expression  of  electric  power  was  deduced  in  case  of  the  large-scale  thermoelectric  device  that  consists  of  the  multiple 
cylindrical  tubes  like  roll  cake  exposed  to  two  thermal  fluids.  The  output  powers  of  the  four  arrangements  were  mathematically  solved  from 
heat  transfer  theory.  The  maximum  output  power  of  roll-cake  module  was  practically  the  largest  for  the  system  of  the  counter  flow  with  a 
single  thermoelectric  tube  (VIC  system),  but  1/4  of  the  ideal  output  power  (VIII  system).  The  multiplication  of  thermoelectric  tubes  can 
shorten  significantly  the  device  length,  although  the  output  power  decreases  only  a  little.  For  example,  the  double  vortical  tubes  (V2CC 
system)  can  generate  95.7%  output  power  by  38.5%  diameter,  compared  with  the  single  tube  (VIC  system). 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  thermal  energy  is  supplied  as  a  form  of  fluid  in  a 
large-scale  thermoelectric  power  plant.  When  a  hot  fluid  of¬ 
fers  the  heat  to  the  hot  junctions  in  the  thermoelectric  mod¬ 
ules,  the  thermoelectric  motive  force  (EMF)  due  to  Seebeck 
effect  is  generated  depending  on  a  temperature  difference 
between  the  hot  and  cold  junctions.  EMF  in  the  open  circuit 
is  the  sum  of  multiplication  of  the  relative  Seebeck  coeffi¬ 
cient  and  the  temperature  difference  Ad  over  all  the  serial 
connections.  In  order  to  gain  the  larger  power,  therefore,  we 
give  the  larger  A 9  to  all  the  thermoelectric  panels  existing 
between  two  hot  and  cold  fluids  [1-4]. 

The  heat  applied  from  a  hot  fluid  is  transferred  firstly  to 
the  panel  surface,  then  into  the  thermoelectric  panel,  and 
finally  to  the  cold  fluid  at  another  surface.  The  two  fluids 
are  resultantly  warmed  or  cooled  along  the  flow  paths,  and 
their  temperature  profiles  through  the  path,  T{x),  change  as 
a  function  of  position,  x.  This  problem  is  known  in  heat 
exchanger  through  an  isolator  [5-7]. 

The  previous  studies  showed  that  the  maximum  of  output 
power  exists  corresponding  to  a  certain  module  size  [8-12]. 
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It  is  because  the  longer  serial  connections  of  thermoele¬ 
ments  can  generate  the  higher  voltage,  but  because  the  out¬ 
put  power  decreases  both  by  the  higher  electrical  resistance, 
and  by  the  smaller  temperature  difference  AT  due  to  the 
higher  heat  transfer  through  the  modules. 

Our  previous  studies  showed  that  the  necessary  panel  area 
could  be  minimized  when  the  thermoelectric  multi-panels 
were  connected  three-dimensionally  [11,12].  The  helical 
system  shown  in  Fig.  1  enables  to  shorten  the  panel  length 
and  makes  the  system  more  compact  and  suitable  for  the 
large-scale  power  generation.  For  example,  the  system  of 
Fig.  1  consists  of  five  thermoelectric  panels,  among  which 
two  hot  fluids  and  the  two  cold  fluids  flow  as  counter  flow. 
The  cylindrical  tubes  consisting  of  the  thermoelectric  mod¬ 
ules  are  suitable  for  circulation  of  the  industrial  fluids,  such 
as  exhaust  hot  gas  or  coolant  water.  Fig.  2  was  the  idea 
that  made  it  practical  to  circulate  the  fluids  as  proposed  in 
Fig.  1  [12].  We  reported  the  mathematical  expression  of  its 
output  power  using  the  heat  transfer  analysis.  However,  the 
practical  arrangement  of  the  multi-tubes  causes  the  practical 
problems  at  the  parts  of  tube  connection,  as  shown  in  Fig.  2. 

Fig.  3  shows  another  possibility  to  shorten  the  cylindrical 
thermoelectric  tubes.  The  two  fluids  flow  in  the  counter  di¬ 
rections  along  the  thermoelectric  vortical  walls,  which  looks 
like  roll  cake.  The  purpose  of  this  work  is  to  deduce  the 
mathematical  expression  of  the  electric  power  extractible 
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Fig.  1 .  Helical  system  with  flat  thermoelectric  panels. 
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Fig.  2.  Cylindrical  multi-tube  system  with  helical  paths  in  which  two 
fluids  flow  in  the  counter  directions  (T3CH  type). 


from  the  thermoelectric  cylindrical  tubes  in  case  of  Fig.  3 
(Vortical  system,  V). 


2.  Basic  assumptions  and  models 

2.1.  Directions  and  paths  of  fluids 

The  previous  studies  clarified  that  the  parallel  flow  did 
not  give  the  larger  output  than  the  counter  flow,  either  in 
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Fig.  3.  Example  of  three-dimensional  helical  system  (roll  pan-cake  model). 
V2CC  system. 


the  flat  panel  systems  (F  system)  [11]  nor  the  cylindrical 
tube  systems  (T  system)  [12].  Therefore,  the  counter  flow  is 
mainly  studied  for  the  vortical  systems  (V  system). 

For  easy  mathematical  handling,  the  roll  cake  systems 
(Figs.  4(a)  and  5(a))  are  simplified  to  the  systems  with  the 
combined  multi-cylinders  (Figs.  4(b),  (c)  and  5(b)).  When 
the  inner  radius  of  the  inner  tube,  r\  \ ,  is  enough  large  com¬ 
paring  with  the  path  width,  a,  we  may  neglect  the  turbu¬ 
lent  effects  near  the  inlet  and  outlet.  The  analyzed  vortical 
(V)  systems  are  named  after  the  number  of  thermoelectric 
tubes  (n)  and  the  direction  of  fluid  flow  (counter  flow,  “C”, 
and  parallel  flow,  “P”).  Therefore,  we  call  the  three  systems 
(Figs.  4(b),  (c)  and  5(b))  as  VIC,  VIP  and  V2CC  system, 
respectively.  For  comparison,  the  case  that  the  both  surfaces 
of  a  tube  are  kept  isothermally  (“I”)  is  also  calculated  as 
VIII  system.  The  hydrodynamic  conditions  and  the  bound¬ 
ary  conditions  in  these  four  mathematical  models  are  listed 
in  Table  1. 
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Fig.  4.  Vortical  systems  with  a  cylindrical  thermoelectric  tube:  (a)  roll-cake  model,  (b)  analyzed  module  (VIC),  (c)  analyzed  module  (VIP). 
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Fig.  5.  Vortical  systems  with  double  cylindrical  thermoelectric  tubes:  (a)  roll-cake  model,  (b)  analyzed  module  (V2CC). 


Table  1 

Conditions  for  the  analyzed  systems 


System 

Hydrodynamic  conditions 

Boundary  conditions 

Model 

VIC 

M\  >  0,  M2  <  0 

7i(0)  =  Tf,  r2( 2tt)  =  Tf 

Fig.  4(b) 

VIP 

Mi  >  0,  M2  >  0 

7i(0)  =  Tf,  72(0)  =  7)" 

Fig.  4(c) 

VIII 

M\  =  00,  M2  =  —00 

Ti(<p)  =  Tf,  T2(<p)  =  7“  at  any  ip 

V2CC 

M\  >  0,  M2  <  0,  M3  =  Mi  >0 

7i(0)  =  7)",  T2(2tt)  =  If,  73(0)  =  T\ (2tt) 

Fig.  5(b) 

2.2.  Cylindrical  thermoelectric  tube 

It  is  assumed  that  all  of  our  thermoelectric  modules  are 
cylindrically  fan-shaped,  and  that  they  consist  of  the  ther¬ 
moelements  with  a  single  layer  of  FI-type  p-n  junctions,  as 
illustrated  in  Fig.  6.  The  thermoelements  are  homogeneously 
aligned  parallel  to  heat  flow,  combined  tightly  without  open 
space,  and  connected  electrically  in  series.  The  hot  and  cold 
fluids  are  isolated  by  this  thermoelectric  panel,  and  they  flow 
along  the  both  cylindrical  surfaces  of  the  panel  in  the  di¬ 
rection  of  cp.  Practically,  as  shown  in  Fig.  6(a),  there  exist 


the  electrodes  connecting  the  thermoelements,  the  insulator 
between  the  elements,  the  protective  insulating  film  on  the 
electrodes,  and  the  fins  that  enhance  the  heat  exchange.  Here 
for  simplicity,  we  neglected  these  inhomogeneous  effects  on 
heat  transfer,  and  considered  them  only  as  the  heat  transfer 
coefficient,  li,  of  the  homogeneously  smooth  modules. 

The  p-n  pairs  of  length,  L,  are  aligned  as  illustrated  in 
Fig.  6(b),  in  order  to  consider  the  heat  transfer  through  the 
thermoelectric  panel  that  depends  on  the  angle,  cp.  This  def¬ 
inition  of  alignment  is  deferent  with  that  of  the  previous 
approach  [12]. 


(a)  (b) 

Fig.  6.  Modeling  of  thermoelectric  modules:  (a)  practical  modules,  (b)  analyzed  modules. 
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3.  Equations  for  output  power 

3.1.  Deduction  of  equations 

Our  previous  analysis  for  multi-flat  panels  [11]  is  fun¬ 
damentally  applied  for  the  multi-tubes.  First  we  examine 
the  heat  transfer  through  the  cylindrical  thermoelectric  tube 
[10,12].  The  heat  transfer  coefficient  K,  is  calculated  by 
building  the  local  heat  balance  equations,  considering  that 
the  heat  transfer  is  not  uniform,  i.e.  the  fluid  temperature 
depends  on  angle,  cp.  Kj  is  deduced  here  as  the  over-all  heat 
transfer  coefficient  per  unit  length  through  the  tube  i  in  the 
direction  perpendicular  to  r  axis, 

1 

( 1  /hp i rp i)  +  ( In (jj  i+i/rUi)/X)  +  (l//q, i+ 1 ru i+  x ) ' 

(1) 

r,  j  and  /;,  y  are  the  radius  of  the  tube  and  the  heat  transfer 
coefficient  between  the  tube  i  and  the  fluid  respectively, 
and  A  is  the  average  heat  conductivity  of  the  module,  defined 
by 

Ap  x  Lp  -F  An  x  Ln  Ap  x  Lp  -f-  An  x  Ln 

A  =  - =  - ,  (2) 

Lp  +  Ln  L 

where  Ap,  An,  Lp  and  Ln  are  the  heat  conductivity  and  the 
length  of  p-  and  n-type  elements,  respectively.  L  is  the  unit 
length  of  thermoelectric  pair,  as  shown  in  Fig.  6(b).  There¬ 
fore,  L  =  Lp  +  Ln.  Note  that  K  has  been  normally  defined 
as  the  over-all  heat  transfer  coefficient  per  unit  length  of  the 
tube  when  the  fluid  flows  along  the  tube  length.  K ,  in  our 
case  is  deduced  when  the  fluid  flows  around  the  tube.  As 
shown  after  the  calculations,  however,  K,  can  be  replaced 
by  KUtt  using  the  normal  formula  K. 

Second,  when  the  system  consists  of  n  tubes,  a  set  of  n  + 1 
simultaneous  derivative  equations  can  be  written  from  the 
heat  transfers  through  tubes  to  obtain  the  fluid  temperature, 
T(cp),  as  a  function  of  angle,  <p.  The  homogeneous  tempera¬ 
ture  along  the  tube  length  is  assumed  that  the  turbulent  flow 
is  well  developed  inside  the  long  path. 


although  this  overestimation  is  not  significant,  as  we  will 
report  separately. 

Third,  the  temperature  dij((p)  at  the  surface  of  the  panel 
i  facing  to  the  fluid  j  is  then  written  using  the  solutions  of 
simultaneous  derivative  equations,  7)(<p)  for  the  angle  (p 

Ou+fw)  ~  Oi.i(cp)  =  y[7i+ 1(^)  -  Tf  cp)}  In  j  .  (6) 

The  electromotive  force,  E ,  is  the  summation  of  the  temper¬ 
ature  difference  over  all  the  panels. 


panel 

£=  E 


i=i 


NfpHxWa 

2: r 


(QiA<p)  -0;+u(<p))d<p 


(7) 


where  a  is  the  difference  between  the  Seebeck  coefficients  of 
p  and  n  elements,  A);,  and  nx  are  the  number  of  thermoelectric 
pairs  in  a  cylindrical  circulation,  and  the  number  densities  of 
pairs  in  the  direction  x,  respectively.  Therefore,  Nq  =  2 n/4> 
and  nx  =  1  /L.  <p  is  the  circumferential  angle  of  a  pair,  as 
shown  in  Fig.  6(b). 

The  electric  resistance,  R,,  is  deduced  by  integrating  the 
resistance  of  a  thin  fan-shaped  layer  in  the  direction  of  ra¬ 
dius. 


Nf-nxw  /  pp  pn  \  |n  / 1  \ 

2tc  \Lp  Ln)  \  rj  j  ) 


(8) 


where  pp  and  pn  are  the  specific  resistivity  of  p  and  n  ele¬ 
ments.  Eq.  (8)  neglects  the  resistance  of  electrode  connect¬ 
ing  p  and  n  elements. 

Finally,  the  output  power,  P .  is  optimized  by  balancing 
the  internal  and  external  resistance. 


(9) 


This  work  uses  hereafter  this  P.  where  the  internal  resistance 
is  equal  to  the  external  resistance. 


3.2.  Output  power  of  VI C,  VIP  and  VIII  systems 


M[Cpp 


d  Tj(<p) 
d(p 


=  TKiw(.Tl(<p)-T2((p)), 


(3) 


M,C  =  ±Ki-\w{Ti-\{(p)  -  Ti(<p))  T  KMTftp) 

d  cp 

—  Ti+\((p)),  (1  <i<n)  (4) 

dTn+i((p) 

M„+iCp,„+i  "+1  y  -  ±Knw(Tn(<p)  -  Tn+l(<p)),  (5) 

d(p 

where  M,  and  Cp,/  are  mass  flow  rate  and  heat  capacity  for 
fluid  i,  respectively.  The  choice  of  sign  ±  depends  on  the 
path  condition.  For  simplicity,  the  heat  transfer  by  Peltier 
effect  and  the  exothermal  heat  by  Joule  effect  are  neglected. 
This  assumption  leads  to  an  overestimation  for  output  power. 


The  simultaneous  equations  are  solved  first  at  the  VIC 
system,  where  one  sheet  of  panel  is  exposed  in  counter 
flow.  EMF  and  the  output  power  are  generally  calculated  by 
Eqs.  (7)  and  (9)  as, 


£Vic 


nxN^MlM2Cp,lCp,2a(T;f  -  If) ln(n2/nt) 


2jtX 

exp(2itKiwD )  —  1 
MiCp.i  +  M2  Cp,  2  exp  (2jtK\  wD) 


(10) 


Pvic  = 


nx{M[M2Cp'  \  Cp2a{T™  -  T™)}2\n(rn/rn) 


87rwA2((pn/Ln)  +  ( pp/Lp )) 
exp(2itK[wD)  —  1  1  " 

MiCp,i  +  M2Cpp  exp(2.7rKiwD)  | 


(11) 
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where 


D 


1 


1 


(12) 


MiCp.i  M2Cpp 

When  MiCp.i  =  —M2CPi2  =  MCp(> 0),  the  special  care 
was  needed  in  solving  the  differential  equations.  EMF  and 
the  output  power  of  this  case  are. 


£Vic  = 


Pm  ic  = 


n x N (f,M CpKiwaiT™  -  T'f )  ln(r12/ni) 
X(MCP  +  2nK\w) 


Ttnxw{MCvK\a{V"  -  rcm)}2ln(ri2//-n) 


(13) 


(14) 


2  X2((pn/La)  +  ( pp/Lp))(MCp  +  2ttK\w)2  ' 

They  are  equal  to  the  infinite  expressions  of  Eqs.  (11)  and 
(12),  respectively,  whenMiCp.i  approaches  to  —M2Cp2. 

The  simultaneous  equations  are  solved  at  the  VIP  system, 
where  one  sheet  of  panel  is  exposed  in  parallel  flow.  EMF 
and  the  output  power  are  calculated  as 


3.3.  Output  power  of  V2CC  system 

In  the  V2CC  system,  we  set  MiCp.i  =  \M2CPt2\  =  MCp. 

£V2CC  =  ~ N(j,nxMCpa(Tf'  -  T “)  sinh(X) 

■  (KiKai  +  K2Ka2)cosh(X) 

+  fKMKg i  -  Kal)  sinh(X) 

yfKfKf  cosh(2Z)  +  (K]  +  K2)  sinh(2X) 


(18) 


Pm 2cc  = 


nxLnLp{MCva(Tla  -  T f )  sinh(X)}2 
2Xw(Lppn  -f-  Lnpp)(Kg\  T  Ka 2) 

■  (KiKal  +  K2Ka2)cosh(X) 

+VKfT2(Ka{  -  Ka 2)sinh(X) 

sfKfKi cosh(2X)  +  (K\  +  /f2)  sinh(2X) 


(19) 


Evip  = 


nxN^M{M2CppCp,2aKfT™  -  Tf)\n(rn/rn) 


2jtX 

1  —  exp(— 2ttK[  wD)  I 

MiCp.i  +  M2CP'2  j  ’ 


(15) 


Pvic  = 


nx{M[M2CpACppMT/'  -  T/WMm/rn) 


1 


SitwX2((pn/Ln)  +  (Pp/Lp)) 
exp(—2TtKiwD)  | 2 


(16) 


MiCp.i  +  M2Cp<2  \ 

Resultantly,  both  the  output  power  /%■  1  c  and  Pvip  are 
independent  from  the  number  of  thermoelectric  elements, 
N<j).  Second,  the  above-mentioned  results  are  deduced  for 
the  situation  that  the  hot  and  cold  fluids  flow  inside  and 
outside  the  thermoelectric  tube,  respectively,  as  illustrated  in 
Fig.  4(b)  and  (c).  However,  these  expressions,  £Vic-  Pm ic> 
Tv  1  p  and  Pyip,  are  the  same  for  the  opposite  situations  that 
the  cold  and  hot  fluids  flow  inside  and  outside  the  tube, 
respectively.  Only  a  difference  is  the  sign  of  EMF. 

When  both  fluids  are  heated  and  cooled  by  the  infinitely 
large  heat  sources,  Tfcp)  and  T2((p)  are  constantly  equal  to 
r}n  and  Tf ,  respectively,  at  any  portion  of  the  tube.  This 
means  that  the  flow  rates  of  both  fluids  become  infinite,  i.e. 
M\  —  \M2\  =  00.  Mathematically,  Pym  is  calculated  from 
Eq.  (9). 


Pvm  = 


Ttnxw{K{oi(T]f  -  Tcln)}2ln(ri2//-ii) 
2  X2((pn/Ln)  +  (Pp/Lp))  ' 


(17) 


where  Ka  1,  Ka2  and  X  were  defined  as 


Kgi  —  — -  In 

71 A 


Kai  =  — ~  In 
ttX 


71  W^/K[K2 
MC~P 


(20) 

(21) 


respectively.  The  definitions  of  Ka\  and  Ka2  are  identical 
with  the  previous  work  [12].  The  output  power  Py2cc  is 
also  independent  from  the  number  of  thermoelectric  ele¬ 
ments,  N(j).  The  obtained  expressions  of  the  output  power 
contain  a  common  term  such  as  nxLnLp{MCpa(T ™  — 
r-)}2/X{LpPn  +  Lnpp ),  however,  the  other  parts  depend 
complicatedly  on  the  geometry  of  the  system,  electric  resis¬ 
tivity  and  heat  transfer  coefficient  K,.  Especially  K,  depends 
on  the  radius  of  thermoelectric  panel  in  the  cylindrical 
system.  This  is  unable  to  simplify  our  expressions,  while 
the  output  power  of  the  flat  multi-panels  systems  could  be 
expressed  by  the  non-dimensional  parameters  [11].  There¬ 
fore,  the  obtained  output  power  for  four  systems  will  be 
numerically  compared  using  the  thermophysical  values  of 
Bi2Te3  semiconductors. 


4.  Physical  properties  and  conditions 

The  thermophysical  properties  of  the  BiTe  thermoelectric 
elements  in  the  literature  scattered  due  to  the  impurities  and 
the  difference  in  preparation  [1].  Table  2  shows  a  typical  set 


Table  2 

Specific  values  of  thermoelectric  elements  at  room  temperature  [1] 


Materials 

Seebeck  coefficient 

Resistivity 

Thermal  conductivity 

Figure  of  merit 

a  (p-V/K) 

p  (|x£2m) 

X  (W/Km) 

Z  (1/K) 

Bi-54.3  at.%  Te  (p) 

162 

5.55 

2.06 

2.605  x  10“3 

Bi-64.5  at.%  Te  (n) 

-240 

10.1 

2.02 
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Table  3 

Parameters  for  thermoelectric  power  generation  system 


Variables 

Values  used  in  this  work 

Thermoelectric  device 

Length,  w 

Thickness  of  tube,  d 

Inner  radius  of  inner  tube,  r\\ 

Distance  between  two  tubes,  a  (thickness  of  fluid  path) 

Number  density  of  pairs,  nx 

Number  of  pairs  round  a  circumferential  cycle,  Nq 

Variable  (typically  2  m) 

0.05  m 

Variable  (typically  2.347  m) 
0.1m 

50  m-1 

100 

Thermal  sources 

Hot  source 

Cold  source 

N2  gas  (inlet  Tj"  =  500  K) 

N2  gas  (inlet  7)"  =  300  K) 

Gas  properties 

Specific  heat,  Cp(hot  =  Cp^ool 

Mass  flow  rate,  Mhot  =  Afcool 

Heat  transfer  coefficient,  /zhot  =  hcoo\ 

1044  J/kgK  (at  400  K,  0.1  Pa) 
Variable  (typically  1  kg/s) 

500  W/m2  K 

of  properties  for  an  identical  sample  at  room  temperature  [  1  ] , 
and  used  for  this  work.  Table  3  shows  the  fluid  properties 
and  the  parameters  for  thermoelectric  tubes.  These  values  are 
identical  with  our  previous  studies  [10-12].  The  thickness 
of  thermoelectric  elements  was  set  as  0.05  m  based  on  the 
previous  discussion  [12].  The  path  width  (the  gap  between 
cylindrical  two  tubes)  was  set  to  be  0.1  m. 


5.  Numerical  comparison  of  systems 

5.1.  Maximum  output  power  in  VIC  system 


The  strict  maximum  of  Eqs.  (11),  (14),  (16)  and  (18) 
are  deduced  as  a  function  of  the  ratio  Ln/Lp,  but  not  taken 
here  because  all  the  equations  hereafter  become  the  longer 
and  complex  mathematical  expressions,  even  if  we  can  ne¬ 
glect  the  dependency  of  K\.  The  output  Pvic  can  be  ap¬ 
proximately  maximized  using  the  well-known  conditions 
that 


X  Pp  -|-  y'V, p  X  pn 

_  s/^n  x  Pp 

^/kn  X  Pp  -|-  yj /. p  X  pn 


and 


(22) 


The  strict  solution  for  the  maximum  was  not  far  from  that 
by  Eq.  (22),  when  we  used  the  numerical  data  mentioned  in 
Section  4.  The  difference  in  output  power  was  only  within 
0.003%,  mainly  because  the  difference  between  kn  and  kp 
was  small.  The  optimized  conditions  for  Ln  and  Lp  (Eq.  (22)) 
are  used  in  the  following  calculations  as  a  good  approxima¬ 
tion. 

Figs.  7  and  8  show  the  output  power  thus  evaluated  for 
VIC  system,  where  the  flow  rate  of  the  fluids,  M\  and  M2,  are 
set  to  be  commonly  1  kg/s  under  the  restriction  of  Eq.  (22). 
The  output  power  has  the  maximum  at  a  certain  length  or 
a  certain  radius.  The  optimum  tube  length  is  given  analyt¬ 
ically  as  wopt-  =  MCp/lnK \  when  the  inner  radius  r\  \  is 
fixed.  For  example.  Fig.  7  shows  that  the  maximum  output 


Fig.  7.  Output  power  from  VIC  system  as  a  function  of  tube  length,  w. 
r\\  is  the  inner  radius. 

power,  Pmax,  is  5.467  kW  at  the  optimum  tube  length  wopt-  = 
2.510m  when  r\ \  =  2m.  When  u)opt-  =  MCp/2ttK\  in  the 
VIC  system  (MiCp,i  =  |M2CPi2|  =  MCp),  Eq.  (14)  gives 
the  maximum  output  power,  pmax. 

pmax  =  nxMCpK{{a(T^  -T^Mm/rn) 

v,c  16k2((pn/Ln)  +  (Pp/Lp)) 

The  optimum  radius  r°p/  when  w  is  fixed,  however,  cannot 
be  solved  analytically.  Numerically  Pmax  =  5.847 kW  is 
obtained  at  /-jj1'  =  2.347 m  when  w  =  2m  (Fig.  8). 


Fig.  8.  Output  power  from  VIC  system  as  a  function  of  inner  radius,  r\\. 
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Fig.  9.  Optimum  tube  length.  u/0P*',  for  the  given  inner  radius,  rn,  and 
mass  flow  rates,  M\  and  M2.  Maximum  output  power  is  shown  only  for 
M  =  1  kg/s.  Open  circles  show  the  optimum  points  shown  in  Fig.  7. 

The  existence  of  the  maximum  output  power  was  also 
reported  in  the  flat  panels  and  the  cylindrical  tubes  [8-12]. 
It  is  qualitatively  because  the  longer  thermoelectric  panel 
can  generate  the  higher  voltage,  but  because  it  increases  the 
electrical  resistance  and  decreases  the  temperature  difference 
A T  between  the  fluids.  It  is  interesting  that  the  values  of 
pmax  are  llearjy  constant  as  shown  in  Fig.  9,  even  if  the 
length  and  radius  are  varied.  This  is  because  K\  changes 
linearly  at  r\  1  >  1  m. 

5.2.  Optimum  dimension  of  VIC  system 

The  relationship  between  the  optimum  length,  u>opt\  and 
the  optimum  radius,  r°j ',  to  obtain  Pnmx  is  analyzed  in 
Fig.  9.  u/opt-  changes  linearly  with  nt,  although  a  small  de¬ 
viation  from  the  linear  regression  is  found  at  r\  \  <  1  m.  This 
linearity  means  that  there  exists  a  certain  optimum  surface 
area,  50pt'  =  27rr^1'  u/opt-  at  a  given  flow  rate.  For  example, 
Sopt-  =  29.8  m2  and  Fmax  =  5.847 kW  for  M  =  1  kg/s  and 
r\  1  =  10  m.  This  constant  .S-opL  shows  again  that  the  cylindri¬ 
cal  tube  can  be  approximated  as  the  flat  panel  at  nt  >  1  m, 
because  the  previous  work  on  the  flat  panel  reported  that  the 
optimum  surface  size  exists  at  the  given  flow  conditions  [11], 

When  r u  is  smaller,  however,  the  flat  panel  approxima¬ 
tion  is  not  valid  because  K\  given  in  Eq.  (1)  does  not  de¬ 
pend  linearly  on  r\\.  For  example,  wopt-  =  51.48  m  and 
Sopt-  =  22.6m2  for  rn  =  0.07  m,  where  Fmax  =  5.807 kW. 
If  we  compare  this  case  with  the  case  of  rn  =  10  m  as 
mentioned  above,  Pmax  decreases  only  by  0.68%  although 
£°pt.  c|ecreases  by  24.2%.  Thus,  the  cylindrical  tube  with 
r\\  <  1  m  can  generate  the  power  more  effectively  per  an 
unit  surface  than  the  flat  panel,  and  we  can  save  the  number 
of  modules  necessary  for  power  generation  using  the  condi¬ 
tion  /'11  <  1  m,  without  dropping  Pmax. 

5.3.  Temperature  profiles  in  VIC  system 

Fig.  10  shows  the  temperatures,  T{(p)  and  0{tp ),  on  the  tube 
surface  in  the  VIC  system.  The  temperature  of  the  fluids 


Fig.  10.  Temperature  profile  along  the  tube  in  VIC  system.  T\  and  T2 
are  temperatures  of  the  fluids,  and  6\  and  62  are  surface  temperatures  of 
the  module,  respectively.  The  inner  radius  and  the  tube  length  are  set  as 
rn  =  2.347  m  and  w  =  2  m,  respectively,  for  the  both  cases  A  and  B. 


Ti(cp)  and  Tiif)  change  linearly  from  the  inlet  to  the  out¬ 
let,  when  MiCp.i  =  lA^Cp^l  =  MCp.  When  the  maximum 
output  is  extracted  from  the  optimum  arrangements  for  the 
length  and  the  inner  radius  (case  A),  the  temperature  dif¬ 
ference  between  T\(q>)  and  Ti{tp )  is  constant,  i.e.  A T((p)  = 
(T'f  —  Tf  )/2,  independent  from  cp,  MCp  and  K\ ,  and  the 
outlet  temperatures  are  T”111  =  =  (Tf  +  Tf1)/!. 

Because  of  the  thermal  resistance  at  the  surfaces,  the  tem¬ 
perature  difference  applied  to  the  thermoelectric  modules, 
A 6((p)  (=  di(ip)  —  02((p)),  is  smaller  than  A T{(p).  As  seen 
in  Eq.  (7),  the  output  power  is  proportional  to  the  integral 
of  A 0{tp),  which  is  visualized  as  the  area  A  surrounded  by 
0\{tp)  and  Oiif)  in  Fig.  10.  When  the  feeding  rate  of  cold 
fluid  I /V/ 2 1  becomes,  for  example,  five  times  larger  (case  B), 
this  area  A  expands  as  shown  in  Fig.  10.  This  is  because 
the  cold  fluid  decreases  the  temperature  02(<Z>).  while  the  hot 
fluid  is  cooled  and  d\((p)  decreases.  In  total,  the  output  power 
increases  when  the  cold  flow  rate  increases. 

Fig.  1 1  shows  the  dependency  of  fluid  flow  rates  on  the 
output  power,  using  Eq.  (11).  Pyic  increases  when  M\  and 


Fig.  11.  Output  power  from  VIC  system  as  a  function  of  flow  rate  of 
hot  fluid,  Mi,  where  ui  =  2m  and  r\\  =  2.347m.  Open  circles  show  the 
data  used  in  Fig.  10. 
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Fig.  12.  Comparison  of  three  systems,  where  the  tube  length  i«  =  2m  and 
the  flow  rate  \M\  =  1  kg/s  commonly.  The  triangles  show  the  positions 
of  maximum  output  power. 


IMtI  increase.  The  largest  Pmax  =  23.392  kW  is  obtained  at 
the  VIII  system  (w  =  2m,  ni  =  2.347 m)  where  the  flow 
rates  of  both  fluids  become  infinite,  i.e.  M\  =  |M2|  =  oo, 
and  T\{cp)  and  Ti(4)  are  equal  to  7’"'  and  T™,  respectively, 
at  any  portion  of  the  tube.  Pvm  in  Eq.  (17)  is  exactly 
four  times  larger  than  P^c  given  in  Eq.  (23)  at  u/>pt' . 
It  is  obviously  because  the  area  A  becomes  two  times 
larger. 

5.4.  Maximum  output  power  in  VIP  and  V2CC  systems 

Fig.  12  shows  the  output  power  for  three  systems,  where 
the  flow  rates  of  the  fluids,  IM/I,  are  set  to  be  commonly 
1  kg/s  and  w  =  2  m.  The  output  power  has  the  maximum  also 
for  VIP  and  V2CC  systems.  Pmax  was  the  largest  for  VIC 
system,  while  the  optimum  inner  radius,  are  1.465  and 
0.903m  for  VIP  and  V2CC  systems,  respectively,  shorter 
than  for  VIC  system  (2.347  m).  The  fact  that  P™q  was 
larger  than  Pyip 's  consistent  with  the  flat  multi-panels  and 
multi -tubes  [11,12], 

Py 2qC  was  95.7%  of  PyiC’  whiie  riiv2cc  was  on^ 
38.5%  of  r°fvlC-  i  inner  surface  of  V2CC  system  per  unit 
tube  length  when  optimized  (2;r(r°1j>t'  +  r)//')  =  12.29  m) 
is  smaller  than  that  of  VIC  (27tr°f  =  14.75  m).  From  the 
viewpoint  of  the  more  compact  system  to  reduce  the  re¬ 
quired  modules  and  the  space,  the  multiple  tube  system  is 
more  effective,  and  the  V2CC  system  seems  the  most  eco¬ 
nomic  for  construction.  This  tendency  is  similar  with  the 
systems  with  the  multiple  flat-panels  or  multi-tubes  [11,12]. 
The  multi-tubes  systems  such  as  V3CC  or  V4CC  systems 
seem  worth  for  further  consideration. 

5.5.  Temperature  profiles  in  V2CC  system 

Fig.  13  shows  the  temperatures,  T(q>)  and  6(4),  on  the 
tube  surface  for  the  case  of  the  maximum  output  in  the 
V2CC  system.  The  temperatures  of  the  fluids  T(tp)  change 


Fig.  13.  Temperature  profile  along  the  tube  in  V2CC  system.  T\,  7?  and 
Tt,  are  temperatures  of  the  fluids,  and  0\\,  0 12,  On,  and  #23  are  the  surface 
temperatures  of  the  module.  ui  =  2m  and  r\\  =  2.347 m. 

exponentially  from  the  inlet  to  the  outlet,  even  if  M\  Cp  i  = 
|M2Cp,2|  =  M^Cpp,  =  MCp.  The  hot  fluid  2  is  cooled 
from  the  both  sides  by  the  coolant,  but  the  outlet  temper¬ 
ature  7 2(0)  for  V2CC  is  higher  than  that  for  VIC  (400 K, 
see  Fig.  10).  Because  the  preheated  coolant  in  the  path  1 
is  again  circulated  to  the  path  3,  the  temperature  increment 
in  the  path  3,  —  73(C)),  is  not  so  significant  as  that 

in  the  path  1,  T\  (2jr)  —  T\  (0).  The  area  Ai  surrounded  by 
012 (cp)  and  611(4)  is  smaller  than  the  area  A2  surrounded 
by  622(41)  and  623(4),  although  Ai  +  A2  is  proportional 
to  EMF. 

In  this  study,  we  have  assumed  that  the  whole  panels  are 
electrically  connected  in  series.  However,  we  may  maximize 
additionally  the  output  power  by  connecting  them  separately 
to  the  external  circuits.  This  possibility  will  be  reported  sep¬ 
arately. 


6.  Conclusion 

This  work  studied  the  thermoelectric  power  generation 
with  cylindrical  multi-tubes  in  which  thermoelectric  ele¬ 
ments  were  embedded.  The  fluid  path  connections  were 
planed  to  be  like  roll-cake.  Four  systems  were  analyzed  us¬ 
ing  the  heat  transfer  theory,  and  their  output  powers  were 
expressed  in  the  analytical  mathematic  equations. 

The  maximum  output  power,  Pmax,  was  obtained  at  the 
VIII  system  where  the  both  surfaces  of  the  tube  were  kept 
isothermally  by  the  infinitively  large  heat  sink.  P™q  was 
the  largest  in  the  other  three  realistic  systems,  although  the 
tube  surface  area  necessary  for  the  maximum  output  power 
was  larger  than  VIP  and  V2CC  system.  When  the  radius  of 
VIC  system  r\\  becomes  larger,  its  output  approaches  that 
of  the  flat  panel,  where  an  optimum  surface  area  generates 
the  maximum  output  power.  The  smaller  r\\  and  the  longer 
w  generates  the  power  more  efficiently  in  VIC  system.  The 
V2CC  system  can  generate  the  95.7%  power  of  that  of  V  1C 
system,  and  r°lj’ty2CC  is  38.5%  of  r°jylc. 
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